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Emission of SN 1006 produced by accelerated cosmic rays 
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Abstract. The nonlinear kinetic model of cosmic ray (CR) acceleration in supernova remnants (SNRs) is used 
to describe the properties of the remnant of SN 1006. It is shown, that the theory fits the existing data in a 
satisfactory way within a set of parameters which is consistent with the idea that SN 1006 is a typical source of 
Galactic CR nucleons, although not necessarily of CR electrons. The adjusted parameters are those that are not 
very well determined by present theory or not directly amenable to astronomical observations. The calculated 
expansion law and the radio-, X-ray and 7-ray emissions produced by the accelerated CRs in SN 1006 agree quite 
well with the observations. A rather large interior magnetic field Bd ~ 100 fiG is required to give a good fit for 
the radio and X-ray synchrotron emission. In the predicted TeV 7-ray flux from SN 1006, the 7r°-decay 7-rays, 
generated by the nuclear CR component, dominate over the inverse Compton (IC) 7-rays, generated by the CR 
electrons in the cosmic microwave background. The predicted source morphology in high energy 7-rays roughly 
corresponds to that of the synchrotron emission. The predicted integral 7-ray flux _F 7 oc e" 1 extends up to energies 
~ 100 TeV if CR diffusion is as strong as the Bohm limit. Only if the interior magnetic field is much lower in the 
SNR, Bd ~ 10 fiG, then the observed 7-ray emission can be due to the accelerated electron component alone. 
In this case, not plausible physically in our view, the lowest permissible value of the electron to proton ratio is 
rather high, and the maximum individual energy and total energy content of accelerated nucleons so small, that 
SN 1006 can not be considered as a typical source of the nuclear Galactic CRs. 
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1. Introduction 

In the last years significant efforts have been made to ob- 
tain direct observational evidence whether the Galactic 
cosmic rays (CRs) - relativistic nucleons and electrons 
- are indeed generated in supernova remnants (SNRs). 
The expected 7r°-decay 7-ray emission, produced in nearby 
SNRs by the accelerated protons in their collisions with 
thermal gas nuclei, is marginally high enough to be de- 
tectable by the present generation of imagi ng at mospheric 



Cherenko v tele scopes (e.g. Drury etal. 1994 ; Naito & 
Takahara |l994t Berezhko & Volk |l997|) . Positive results 
of such observations would constitute a necessary condi- 
tion for a dominant role of SNRs in the production of the 
Galactic CRs and of their energy spectrum up to the knee 
energy ~ 10 15 eV. 

If SNRs accelerate nucleons at least to 100 TeV/n as it 
is req uired for Galactic CRs (e.g Berezhko & Ksenofontov 
1999] ) , then in smooth extension of the expected hadron- 
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ically g enera ted 7r°-d ecay 7 -ray peak at 67.5 MeV (e.g. 
Stecker 1971 ; Dermer 1986] ) one expects the nucleonic 7- 
ray spectrum to have an almost single power-law form to 
about 1 TeV, before turning over above 10 TeV. At the 
same time the leptonic Inverse Compton (IC) 7-ray spec- 
trum may be quenched at a lower cut-off energy due to 
electron synchrotron losses (as in our model). Then 7-ray 
observations within the range from 10 MeV to 100 TeV 
would be useful in discriminating between the nucleonic 
and leptonic components. 
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Recent observations of nonthermal X-rays and 7-rays 
indicate that at least CR electrons are accelerated in 
SNRs. SN 1006 is one of the SNRs for which there is 
evidence that electrons reach energies of about 100 TeV 
(Koyama et al. |1995| ; Tanimori et al. |1998] , |2001| ). It is 
also one of the three shell type SNRs in which TeV 7-ray 
emission has been detected up to now. However, the inter- 
pretation of these data is not unique. Depending on the 
assumed values for the unknown physical parameters of 
SN 1006 (mainly the value of the magnetic field, the elec- 
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tron to proton ratio, and to some extent also the nucleon 
injection rate), the observed high-energy 7-ray emission 
can be predominantly either inverse Compton radiation 
due to CR electrons scattering on the microwave back- 
groun d rad iation (as predicted by Mastic hiadis ( 1999 ), 
Pohl ( |L996[ ), Mastic hiadis & de Jager ( |1996| ), and Yoshida 
& Yanagita (1997) from the X-ray synchrotron emis- 
sion, and through the interpretation of the X-ray and 
7-ray data by Aharonian (1999), Aharonian & Atoyan 
(1999), Berezhko et al. (1999)), or 7r°-decay emission due 



t o had ronic collisions of CRs wit h gas nuclei (Aharonian 
( I l999|) , Aharonian & Atoya n (|l999|) , Berezhko et al. 
(|l999|) , Berezhko et al. (|200l] )). 

The thermal and nonthermal X-ray emiss ion has re- 
cently been rediscussed by Allen et al. ( 2001 ) who draw 
phcnomcnological conclusions on the 7-ray emission to be 
dominated by IC radiation as well as on the production 
and the characteristics of the nonthermal nucleonic parti- 
cle component which they estimate to have a total energy 
of 10 50 ergs. 

In contrast, our starting point is the overall SNR dy- 
namics and the acceleration theory of the nonthermal 
component. Thus we accept the X-ray results and the 
quantitative distinction between thermal X-ray emission 
and the nonthermal synchrotron components, including 
the constraints on the external density derived from them. 
We calculate the outer shock size, speed and compression 
ratio as well as the nonthermal quantities as solutions of 
the nonlinear, time-dependent kinetic equations for elec- 
trons and protons as functions of space, time and parti- 
cle momentum, rather than phenomenologically assuming 
forms of the electron and proton momentum distribution 
functions. 

In this way we also obtain full morphological informa- 
tion. We determine the unknown parameters like the up- 
stream magnetic field strength and the electron to proton 
ratio, and constrain the nucleon injection rate which can- 
not yet be accuratetely calculated from known theory, by 
comparing with the observations. For this purpose we use 
the selfconsistent kinetic model of diffusive acceleration 



of CRs in SNRs (Berezhko et al. 1996; Berezhko & V61k 



1997) and investigate the 7-ray emission from SN 1006. A 



preliminary version of our results is given in Berezhko et 
al. ( |200lD . 



In contrast to a previous study (Berezhko et al. 1999) 
we restrict ourselves to the so-called Bohm limit for CR 
diffusion near the shock, assuming efficient and strong ex- 
citation of magnetohydrodynamic waves by the accelerat- 
ing particles themselves. This is consistent with the large 
magnetic field strengths which we infer for this remnant. 

Our considerations show that, together with a renor- 
malization due to the lack of spherical symmetry of the 
nucleon injection, the existing SNR data are consistent 
with very efficient acceleration of CR nuclei at the SN 
shock wave which converts a significant fraction of the 
initial SNR energy content into CR energy as required for 
a typical source of the Galactic CR nuclei. The relative 
amount of energetic electrons is rather small, and may 



require additional sources to SNe type la like SN 1006. 
Therefore the observed 7-ray emission of SN 1006 can in- 
deed be of hadronic origin, and we consider this as the 
physically most plausible solution from the point of view 
of acceleration theory. 

Nevertheless, the existing observations do not strongly 
exclude a solution in which nuclear CRs play no important 
role and all nonthermal emissions are of leptonic origin. 
Therefore we analyze whether the existing data can also 
be fitted by an essentially different set of parameters. We 
demonstrate that there is an alternative possibility to fit 
the overall data, albeit with somewhat lower quality. 

It implies a (physically not plausible) much lower in- 
jection rate and acceleration efficiency of protons, a rather 
large electron to proton ratio, and a lower magnetic field 
value compared with the case of efficient CR nucleon ac- 
celeration. If we would assume SN 1006 to be a typical 
representative, in this second case we could not consider 
the SNRs as the source population of the Galactic CRs 
due to the low proton acceleration efficiency. To discrim- 
inate empirically between these rather different scenarios 
from our results 7-ray measurements above 10 TeV are 
required: in this very high energy range we expect a mea- 
surable 7-ray flux only in the case of efficient nucleonic CR 
production. We also confirm the earlier result (Aharonian 



1999; Aharonian fc Atoyan 1999 Berezhko et al. 1999 



Berezhko et al. 2001) that the IC 7-ray emission should 
cover the SNR almost uniformly, whereas the 7r°-decay 
emission should be peaked behind the shock front and 
limited to two polar caps where the interstellar magnetic 
field is parallel to the shock normal. 

2. Model 

Since SN 1006 is a type la supernova (SN) we suggest that 
its evolution takes place in a uniform interstellar medium 
(ISM), not modified by a wind from the progenitor star. 
The general picture is then well-known. 

A SN explosion ultimately ejects a shell of matter with 
total energy E sn and mass M e - y During an initial period 
the shell material has a broad distribution in velocity v. 
The fastest part of these ejecta is d escrib ed by a po wer law 
dM ei /dv oc v 2 ~ k (e.g. Jones et al. |l98l| ; Chevalier |1982| ). 
The interaction of the ejecta with the ISM creates a strong 
shock there which accelerates particles. 



Our nonlinear model (Berezhko et al. 1996; Berezhko 
& Volk 1997 ) is based on a fully time-dependent solu- 
tion of the CR transport equation together with the gas 
dynamic equations in spherical symmetry. It yields at any 
given phase of the SNR evolution the complete spatial dis- 
tribution of gas and accelerated CRs and their spectrum. 
This makes it possible to calculate any kind of emission, 
and its morphology, produced by the accelerated particles. 

The CR diffusion coefficient is taken as the Bohm limit 



k(p) = K(mc)(p/mc), 



(1) 



where approximately k(toc) = mc 2 / {3eB), e and m are 
the particle charge and mass, p denotes the particle mo- 
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mentum, B is the magnetic field strength, and c is the 
speed of light. With this diffusion coefficient we assume 
that the production of scattering waves is so strong that 
limitations due to wave refraction (Malkov et al. 1995) do 
not come in to influence the solution. 

The number of suprathermal protons injected into the 
acceleration process is described by a dimensionless injec- 
tion parameter r\ which is a fixed fraction of the ISM par- 
ticles entering the shock front. For simplicity it is assumed 
that the injected particles have a velocity four times higher 
than the postshock sound speed. Unfortunately there is 
no complete selfconsistent theory of a collisionlcss shock 
transition which can predict the value of the injection rate 
and its dependence on the shock parameters. For the case 
of a purely parallel plane shock hybrid simulations pre- 
dict a quite high ion in jection (e.g. Scholer et al. 1992 ; 
Bennet & Ellison 1995 ) which corresponds to the value 
r\ ~ 10~ 2 of our injection parameter. Such a high injec- 
tion i s con sistent w ith an alytical models (Malkov & Volk 
1995 , 1996| ; Malkov 1998|) and confirmed by measur ement s 
near the Earth's bow shock (Trattner & Scholer 1994 ). 



We note however that in our spherically symmetric model 
these results can only be used with some important mod- 
ification. In reality we must consider the evolution of the 
large scale SN shock which expands into the ISM with 
its magnetic field. In the case of SN 1006, at the cur- 
rent evolutionary phase, the shock has a size of several 
parsecs. On such a scale the unshocked interstellar mag- 
netic field can be considered as uniform since its random 
component is characterized by a much larger main scale of 
about 100 pc. Then our spherical shock is quasi-parallel in 
the polar regions and quasi-perpendicular in the equato- 
rial region. This magnetic field essentially suppresses the 
leakage of suprathermal particles from the downstream 
region back upstrea m whe n the shock is more and more 
oblique (Ellison et al. |l995fc Malkov & Volk |l995| ). Applied 
to the spherical shock in the uniform external magnetic 
field it would mean that only small regions near the poles, 
covering less than 10% of the shock surface, allow a suf- 
ficiently high injection that ultimately leads to the trans- 
formation of an essential part (more than a few percent) 
of the shock energy into CR energy, whereas the main 
part of the shock is an inefficient CR accelerator. If one 
takes into account the Alfven wave excitation due to CR 
streaming (which becomes efficient already at a very low 
injection rate rj ~ 10~ 7 ) the local injection rate has to 
be averaged over the fluctuating magnetic field directions 
and is lower than for the purely parallel case by a factor of 
hundred. Therefore we adopt here the value r\ ~ 10~ 4 for 
the injection parameter (Volk et al. 2002] ). The detailed 
choice 7] = 2 x 10~ 4 in the next section is made in order to 
achieve a nonlinear shock modification that is consistent 
with the observed radio spectral index. 

According to our above estimate a substantial part of 
the shock still efficiently injects and accelerates CRs. This 
fraction is effectively increased relative to the above per- 
centage due to broadening of the injection region by the 
strong wave field as well as by CR diffusion perpendicular 



to the mean magnetic field. In addition, the overall conser- 
vation equations ensure an approximately spherical char- 
acter of the overall dynamics. Therefore, we assume the 
spherically symmetric approach for the nonlinear particle 
acceleration process to be approximately valid in those 
shock regions where injection is efficient. To take the effec- 
tive injection fraction / re < 1 into account, we then need 
to introduce a renormalization factor for the CR acceler- 
ation efficiency and for all the effects which it produces in 
the SNR. A rough estimate performed for the simple case 
of spherical shoc k expe nded in the uniform outer magnetic 
field (Volk et al. |2002|) gives / re = 0.15 to 0.25. 

Note that such a picture is consistent with the ob- 
served structure of SN 1006: the intense radio and X-ray 
emissions come from two bright rims with radia lly ori- 
ented magnetic field (Reynolds & Gilmore 1993 ). They 
may be the polar regions of a quasispherical shock in an 
outer magnetic field. The 7-ray observations may indicate 
a similar asymmetry (Tanimori et al. 1998) which would 
be well explained by a corresponding concentration of the 
accelerated nucleons towards the polar regions. 

We assume that electrons are also injected into the 
acceleration process, still at nonrelativistic energies be- 
low m e c 2 . Since the electron injection mechanism is not 
very well known (e.g. Malkov & Drury 2001) for simplic- 
ity we consider their acceleration starting from the same 
momentum as that of protons. At relativistic energies (for 
protons) these electrons have exactly the same dynamics 
as the protons. Therefore, neglecting synchrotron losses, 
their distribution function at any given time has the form 



Up) = K ep f(p) 



(2) 



with some factor K cp which is about 10 2 



for the average 
CRs in the Galaxy. We take it here as a parameter. 

The electron distribution function f c (p) deviates only 
at sufficiently large momenta from this relation due to 
synchrotron losses, which are taken into account by sup- 
plementing the ordinary diffusive transport equation with 
a loss term: 

— = v«v /c - „v/. + — p— _ (_/. j , (3) 

where the first three terms in the right hand side of this 
equation describe diffusion, convection due to the mass 
velocity w of the gas and adiabatic effects, respectively. 
The synchrotron loss time in the third term is determined 
by the expression (e.g. Berezinskii et al. 1990| ) 



ArlB 2 p 



(4) 



where m c is the electron mass and ro the classical electron 
radius. 

The solution of the dynamic equations at each in- 
stant of time yields the CR spectrum and the spatial 
distributions of CRs and gas. This allows us to calcu- 
late the expected flux F^(e 7 ) of 7-rays from 7r°-decay due 
to hadronic (p-p) collisions of CRs with the gas nuclei. 
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Following the work of Dermer ( 198E ) and its later improve- 
ment by Naito & Takahara (1994) we use here the isobar 
model at the protons kinetic energies ek < 3 GeV and the 
scaling model at tk > 7 GeV with a linear connection be- 
tween 3 and 7 GeV. This model agrees very well with the 
simpler approach, introduced by Drury et al. (1994) (see 
also Berezhko . 



Volk 



1997, 2000 



Berezhko et al. 



1999) at 



high energies e 7 > 0.1 GeV, except in the cutoff region, 
where the scaling model yields an essentially smoother 
turnover of the 7-ray spectrum at somewhat lower ener- 
gies. 

The choice of K ep allows us to determine in addition 
the electron distribution function and to calculate the as- 
sociated emission. The expected synchrotron flux at dis- 
tance d from th e SNR is given by the expression (e.g. 
Berezinskii et al. 199C ) 



Sv = 3 x * 21 1 00 drr2B± r dpp2f ^ p)F fe) (5) 

in erg/(cm 2 s), where 

poo 

F(x) = x K 5/3 (x')dx', 

J X 

K^(x) is the modified Bessel function, v c = 
3eB±p 2 /[47r(m e c) 3 ] and B± is the magnetic field 
component perpendicular to the line of sight. 

Since in the shock region where particle injection and 
acceleration is efficient the upstream magnetic field is as- 
sumed to be almost completely randomized due to intense 
Alfven wave generation. For such a strongly turbulent up- 
stream field, consistent with Bohm diffusion, the post- 
shock field strength B2 is given by by B2 = a Bo, where 
a% = 1/3+ (2/3)cr 2 . We shall approximately equate ub 
with the gas compression ratio a at the shock that We 
suggest also that in such a young SNR like SN 1006 the 
downstream magnetic field £?d in a relatively thin region 
between the shock and the ejecta is approximately uniform 
and we shall take B& = B(r < R s ) = Bi- 

The relativistic electrons produce 7-ray emission due 
to IC scattering of background photons. It is not difficult 
to show that due to the hard spectrum of accelerated elec- 
trons only the 2.7 K cosmic microwave background (CMB) 
is important in the case considered. The expected differ- 
ential flux of IC 7-rays as a function of their energy e 7 can 
be represented in the form: 



4nc 



drr 2 
"'O 



den ph (e) 



x / dpp 2 a(e c ,e~ / ,e)f c (r,p) 



(6) 



in photons/(cm 2 s erg), where (Blumenthal & Gould 1970) 

, 3cr T (m c c 2 ) 2 
o-(e e ,e 7 ,e) = 



4ee 2 



2q\nq+ (1 +2q)(l - q) +0.5 



(r<?) 2 (i-g) 



l + Tq 



(7) 



is the differential cross section for the up-scattering of a 
photon with incident energy e to energy e 7 by the elastic 
collision with an electron of energy e c , 



1 



n-ph 



n 2 (hc) 3 exp(e/k B T) - 1 



(8) 



is the blackbody spectrum of the CMB, h and ks are 
Planck and Boltzmann constant respectively, T = 2.7 K, 
o"t = 6.65 x 10~ 25 cm 2 is the Thomson cross-section, q = 
e 7 /[r(e c — e 7 )], T = 4ee e /(m e c 2 ) 2 , and p m i n is the minimal 
momentum of electron, whose energy e e is determined by 
the condition q = 1. 

Energetic electrons produce also nonthermal brem- 
strahlung (NB) emission, interacting with both ambient 
electrons and protons. We use standard formulae for the 
cross section to calculate the NB 7-ray flux F^ B (e~ ( ) (e.g. 
Baring et al. 1999). As it was already demonstrated by 



these last authors, in young SNRs the primary energetic 
electron number density due to electrons directly acceler- 
ated from the suprathermal postshock distribution by far 
exceeds that of secondaries which were produced via the 
decay of charged pions created in p-p collisions. Hence 
the contribution of secondaries to the NB, IC, and syn- 
chrotron emission can be neglected. 



3. Results and Discussion 

Since SN 1006 is a type la SN we use typical SN la pa- 
rameters in our calculations: ejected mass M e j = 1.4M©, 
k = 7, and a uniform ambient ISM. 

The ISM density is the most relevant parameter, espe- 



cially for the 7-ray production (e.g. Berezhko et al. 1999). 
Here we use the most appropriate value of the ambient 
number density of ISM hydrogen TVh = 0.3 cm -3 as well 
as the distance d = 1.8 kpc, consistent with X-ray and op- 
tical imagery of SN 1006 (Winkler & Long [l99^ , see also 
Allen et al. 2001 and references therein). 

The value To — 10 4 K is used for the ISM temperature. 
Note that SNR and CR dynamics are not sensitive to To, 
because the shock structure is mainly determined by the 
Alfvenic Mach number. 

The other important parameter is the ISM magnetic 
field Bo- Unfortunately there is no direct way to deter- 
mine its value from the observations. Since it influences 
the synchrotron spectrum in an essential way, we use an 
upstream magnetic field value Bo = 20 /xG. It is required 
to provide the observed shape of the synchrotron spectrum 
in the radio and the X-ray bands. 

The gas dynamics problem is characterized by the fol- 
lowing length, time, and velocity scales: 



Ro = (3M cj /4^ ) 1/3 , t = Rq/Vq, V = ^2E sn /M ei , 

which are the sweep- up radius, sweep- up time and mean 
ejecta speed respectively. Here p — lArripNn is the ISM 
mass density; m p is the proton mass. 



Berezhko et al.: Emission of SN 1006 



5 



The shock expansion law during the free expansion 
phase (t < to) is then 



R s cx £S~ 3)/2k Po 1/k t( k - 3 ^ k - 2 > 
(Chevalier 1982| ) which for k = 7 gives 
R s cx {ElJp,flh^. 
In the adiabatic phase (t > to) we have 
Rs cx (£ sn /A)) 1/5 i 2/5 . 



(9) 



(10) 



(11) 



The observed expansion law of SN 1006 (Moffett et al. 



1993J) is R s cx f with fx = 0.48 ± 0.13. We conclude that 
within the observational errors SN 1006 should already be 
in the adiabatic phase. 

The calculations together with the experimental data 
are shown in Fig.|l]-||. For the assumed SN distance and 
ISM density an explosion energy E sn = 3 x 10 51 erg fits the 
obser ved SN R size R s and its expansion rate V s (Moffct 
et al. 



1992) 



The shown uncertainties of R s and V s are 
calculated as the uncertainties of the angular size and ex- 
pansion rate times the distance d = 1.8 kpc. 

According to Fig]j]a SN 1006 is indeed already in the 
adiabatic phase. The assumed injection rate ?) = 2x 10~ 4 
leads to a significant modification of the shock which at 
the current epoch, t — 995 yr, has a total compression 
ratio a = 6.3 and a subshock compression ratio a s = 3.6 
(Fig.|b). 

The different SNR energy components, the ejecta en- 
ergy E e j, the gas kinetic (i?gk), and thermal (E gt ) ener- 
gies, and the CR energy E c , are presented in Fig.|l]c as 
functions of time. The acceleration process is character- 
ized by a high efficiency in spherical symmetry: at the 
current time t/to — 4.68 about 53% of the explosion en- 
ergy has been already transferred to CRs, and the CR 
energy content E c continues to increase to a maximum of 
about 60% in the later Sedov phase (Fig.[j]c), when par- 
ticles start to leave the source. As usually predicted by 
the model, the CR acceleration efficiency is significantly 
higher than required for the average replenishment of the 
Galactic CRs by SNRs, corresponding to E c « 0.1-E sn - As 
discussed before, this discrepancy can be attributed to the 
physical conditions at the shock surface which influence 
the injection efficiency. The magnetic field geometry is the 
most important factor: at the quasiperpendicular portion 
of the shock ion injection (and subsequent acceleration) 
is presumably depressed compared with the quasiparallel 
portion. Therefore the number of nuclear CRs, calculated 
within the spherically-symmetrical approximation, should 
be renormalized by this depression factor (see Volk et al. 



2002 for details). Assuming SN 1006 to be an average 
Galactic CR source, the renormalizing factor should be 
/ re = 0.2. It means that according to our calculation the 
actual energy of accelerated protons in SN 1006 at the 
present epoch is 
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Fig. 1. (a) Shock radius R s and shock speed V s ; (b) total 
shock (a) and subshock (er s ) compression ratios; (c) ejecta 
(i? c j), CR (E c ), gas thermal (E gt ), and gas kinetic (2£ g k) 
energies as a function of time. Scale values are Rq = 3.2 pc, 
Vo — 14675 km/s, t = 212 years. The dotted vertical line 
marks the current epoch. The observed size and speed of 
the shock (Moffett et al. |l993j ), are shown as well. 



Note that such a factor is consistent with observations 
which show that efficient CR acceleration takes place at 
about 25% of the shock surface (e.g. Allen et al. 2001). 
The volume-integrated (or overall) CR spectrum 



(13) 



E c = 0.53/ rc £ sn = 3.2 x 10 au erg 



(12) 



N(p,t) = 16ttV / drr 2 f(r,p,t) 
Jo 

has, for the case of protons, almost a pure power-law form 
N oc p -7 over a wide momentum range from 10 _2 m p c up 
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log(p/m p c) 

Fig. 2. The overall (= spatially integrated) CR spectrum 
as function of momentum. Solid and dashed lines corre- 
spond to protons and electrons, respectively. Thick and 
thin lines correspond to efficient and inefficient proton ac- 
celeration, respectively. 

to the cutoff momentum p max = e max /c, where e max » 
3 x 10 14 eV is the maximum CR energy (Fig.||). The lim- 
itation of this numerically determined value p max can be 
understood mainly in terms of the influence of geometrical 
factors which are the finite size and speed of the shock, its 
deceleration and the adiabatic cooling effect in the down- 
stream region (Berezhko 1996 ). Due to the shock modifi- 
cation the power-law index slowly varies from 7 = 2.1 at 
p ^5 m-pC to 7 = 1.8 at p > 10 3 to p c. 

The shape of the overall electron spectrum N c (p) de- 
viates from that of the proton spectrum N (p) at high mo- 
menta p > p\ w 10 3 rn p c, due to the synchrotron losses in 
the downstream region with magnetic field k, 120 fiG 
which is assumed uniform in this region (B^ rj B 2 — &B ). 
According to expression (Q) the synchrotron losses should 
become important for electron momenta greater than 



Pi 



1.3 



10 8 yr\ /10 
t J { B d 



(14) 



Substituting the SN age t = 10 3 yr into this expression, 
we have p\ « 600m p c, in reasonable agreement with the 
numerical results (FigJ|). 

In detail the momentum dependence of the spatially 
integrated electron spectrum N e from Eq. (13) is some- 
what more complicated at high energies. This can be un- 
derstood as follows: 

The shock constantly produces the electron spectrum 
f c oc p~ q with q w 4 up to the maximum momentum p° nax 
which is much larger than p\. This spectrum is not mod- 
ified by synchrotron losses within the downstream region 
of thickness 



I = T1U2, 



(15) 



where u 2 = V s /a is the gas speed with respect to the 
shock. At low momenta p < p\ the loss time T\ exceeds the 



age of the system. Formally this means that the size of the 
lossfree region is greater than the thickness Ar « i? s /(3cr) 
of the swept up shell. At large momenta p > p\ we have 
I < Ar. Within the downstream region r < R s — I the 
electron spectrum therefore becomes essentially modified 
(f c oc p~ 5 ) by the influence of the synchrotron losses. At 
a given radius r < R s accelerated electrons have an age of 
At « (Rs — r)/u2 and therefore, cf. Eq. (14), their spec- 
trum has a break at p\(r) — p\ x (t/At). According to 
its definition ([bj]) the spatially integrated electron spec- 
trum N c (p) is a mixture of electron spectra with different 
p\ (r) . Qualitatively one would expect that within the mo- 
mentum range p\ to p^ lax , the overall electron spectrum is 
close to the form N e oc p~ 3 , somewhat hardening towards 
Pmax- This is in agreement with the calculation. 

The maximum electron momentum p^ax can be esti- 
mated by equating the synchrotron loss time (^) and the 
acceleration time 



3 

Au 



Hi 



'21 

U 2 



(16) 



where u\ = V s and u 2 = u\ja are the upstream and 
downstream gas velocities relative to the shock front, and 
Am = u\ — U2- During the acceleration time r a the elec- 
trons spent the time 3ki/(miAm) in the upstream mag- 
netic field Bq, and 3k2/(u2Au) in the downstream mag- 
netic field B2 = (JbBq (in the case considered ctb = er). 
Therefore, in the general case the maximum electron mo- 
mentum p^ax is a solution of the equation 



-1 



rf 1 ' 



(B a , p) +rf 1 (B 2 , p)a/a B ]/(l + a/a B ), (17) 



which can be written in the form 

T) C / V 



6.7 x 10 4 



10 3 km/s 



((7-1) /10/xG 



er(l + (TbO") \ Bq 



(18) 



The main part of the electrons with the highest energies 
e <^ 10 TeV is produced at the end of the free expansion 
phase. At this stage V s ~ Vq which leads to a maximum 
electron momentum p^, ax ~ 4 x 10 4 m p c in agreement with 
the numerical results (Fig.||). 

The parameters K cp = 3 x 10~ 4 and B± « 0.3i?d ~ 
36 pG provide good agreement between the calculated and 
the measured synchrotron emission in the radio- and X-ray 
ranges (Fig.||). The steepening of the electron spectrum 
at high energies due to synchrotron losses naturally yields 
a fit to the X-ray data with their soft spectrum. Note 
that two kinds of X-ray data are presented in Fig.||. The 
Hamilton et al. (1986) data represent the total (i.e. ther- 
mal plus nonthermal) X-ray emission in a relatively wide 
frequency range, where the nonthermal synchrotron emis- 
sion contribution dominates at frequencies v <; 1 18 Hz. 
At these highest frequencies the Allen et al. (|l999|) RXTE 
nonthermal data are also presented (in red color), lying 
slightly below the Hamilton et al. data points. Clearly our 
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nonthermal model spectra must lie below the Hamilton et 
al. spectral points at v < 10 18 Hz. Such a smooth spectral 
behavior is achieved in a 20 /iG upstream field. 

Since the total number of accelerated protons has to 
be reduced by the factor / rc = 0.2, the same number of 
accelerated electrons then corresponds to a renormalized 
parameter K cp — 1.5 x 10" 3 . Then the total energy of 
accelerated electrons in SN 1006 is = 5 x 10 47 erg. 

The ratio of accelerated electrons to protons thus 
comes out lower by a factor of about 6 than the canonical 
ratio of 0.01, generally deduced from observations of the 
Galactic CRs. This is an interesting consequence of our 
model: if SN 1006 is typical for the nuclear Galactic CR 
source population, then other sources like young Pulsars 
or Pulsar Nebulae must also significantly contribute to 
the G alactic CR electr on po pulation (e.g. Aharonian et al. 



1995; Pohl & Esposito 1998). An additional fraction might 



come from secondary interstellar electrons/positrons due 
to the decay of charged pions that are produced in inelas- 
tic collisions of CR nuclei with nuclei of the interstellar 
gas. Since the majority of Galactic Supernovae is not of 
the type la to which SN 1006 belongs, but rather occurs 
as a consequence of core collapse events, there is the alter- 
native possibility that core collapse Supernovae produce 
a higher electron to proton ratio. Such an increase might 
have to do with the different circumstellar magnetic field 
structure of massive stars, but this is indeed only a possi- 
bility for which we have no proof at present. 

The radio data are fitted by a power law spectrum 
S v oc v~ a , whose index a = 0.57 ±0.06 (Allen et al. |200ip 
is noticeably larger than 0.5, corresponding to the electron 
spectrum N c oc p~ 2 produced by an unmodified shock 
with compression ratio a = 4. In our case the shock is es- 
sentially modified by the backreaction of the accelerated 
protons (see Fig.[l|b): its compression ratio a = 6.3. At the 
same time low energy electrons with momenta p ^ 10m p c 
(e c ^ 10 GeV), which produce synchrotron emission at 
v < 10 GHz, are accelerated at the subshock which has 
the compression ratio a s = 3.6. Therefore these electrons 
have a steeper spectrum N c oc p~ 21 that leads to the 
expected radio spectrum S v oc i/ -0 - 54 , fitting the experi- 
mental data very well (FigJ|). The fact that the observed 
value of the radio power law index a exceeds the value 
0.5 can be considered as an indication that the shock 
is essentially modified. A relatively high upstream mag- 
netic field strength Bq — 20 //G, compared with typi- 
cal ISM values Bq ~ 5 fj,G, and a corresponding down- 
stream value B±_ = 36 fj,G are required to have radio 
emitting electron energies in the steep part of their spec- 
trum e c < 1 GeV on the one hand, and to give a smooth 
rollover of the synchrotron spectrum S v (v) at frequencies 
v = 10 15 to 10 18 Hz on the other. According to model cal- 
culations by Lucek and Bell ( 20001 ), the existing ISM field 
can indeed be significantly amplified near a strong shock 
by CR streaming. 

Such a strongly nonlinear behavior has three conse- 
quences in that (i) it leads to the high magnetic field 
strength, required to understand the measured nonther- 




Fig. 4. IC (dashed lines), 7r°-decay (solid lines) and NB 
(dash-dotted lines) integral 7-ray energy fluxes as a func- 
tion of 7-ray energy for the same cases as in Fig.^. High 
energy 7-ray data (Tanimori et al. 199S ) , the E GRET up- 
per limit as given by Mastichiadis & deJager 1996 , and 
the JANZOS upper limit (Allen et al. [1995D are shown. 



mal synchrotron spectrum, (ii) it is consistent with our 
assumption of Bohm diffusion, and (iii) it suggests a value 
7] r~ 10~ 4 for the injection parameter that we use (see also 
section 2). 

The calculated integral IC and 7r°-decay 7-ray energy 
fluxes are presented in Fig.|| and || together with the avail- 
able experimental data. 

In Fig. 4 we present also the integral NB 7-ray fluxes; 
they play no important role for SN 1006. Note that the 
number of accelerated protons, which produce 7r°-decay 
7-rays, has been reduced by the factor / re = 0.2. 

According to the calculation, the hadronic 7-ray pro- 
duction exceeds the electron contribution by a factor of 
about 7 at energies e 7 ^5 1 TeV, and dominates at e 7 > 
10 TeV (FigJ||). The calculation is in good agreement with 
the TeV-measurements rep orted by the CANGAROO col- 
laboration (Tanimori et al. 1998 ), and it does not contra- 
dict the EGRET upper limit F E = 8 x 10~ 9 cm^s" 1 at 



e 7 = 1.4 GeV (cf. Mastichiadis & de Jager |996|. This is 
also confirmed by Figj5[ where we compare our calcula- 
tions with the revised CANGAROO data (Tanimori et al. 
|2001| ) plus an extended set of EGRET upper limits (Naito 
et al. |1999| ). 

The differential 7r°-decay 7-ray flux dF^ /de 7 has the 
expected peak at photon energy = 67.5 MeV (see 

also Naito & Takahara 1994). This feature may be used 
for the identification of a hadronic contribution to the ob- 
served 7-ray flux at the particle energies < 10 GeV mainly 
producing it. However, two considerations weaken this 
possibility. First of all, for reasons of instrumental sensi- 
tivity, we expect only nearby SNRs to be detectable in the 
foreseeable future. Such objects have diameters approach- 
ing 1 degree, and SN 1006 is an example. The diffuse 
Galactic 7-ray background is large at these energies due to 
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Fig. 5. Differential 7r°-decay (solid lines) and IC (dashed 
lines) 7-ray energy fluxes as a function of 7-ray energy 
for the same cases as in Fig.||. The recent differ ential high 
energy 7-ray energy flux data (Tanimori et al. 2001 ) and 
EGRET upper limits (Naito et al. 199S) are also shown. 



the steep spectrum of the Gal actic CRs that produce it, as 
pointed out by Drury et al. ( 1994 ). Therefore, we expect 
SNRs to be best observable at much higher 7-ray energies 
e 7 > 100 GeV, where the hard source spectra typically 
dominate over the background emission. Secondly, from a 
strictly empirical point of view, such a nuclear CR com- 
ponent that produces the 7r°-decay bump need not extend 
far beyond the GeV region. For the existence, in a SNR, of 
a nuclear CR component approaching the energy region of 
the knee in the Galactic spectrum, the detection of a 67.5 
MeV 7-ray feature is therefore neither a necessary nor a 
sufficient condition. 

On the other hand, due to the high magnetic field 
strength which we deduce for this remnant, the 7-ray 
spectra produced by the electronic and hadronic CR com- 
ponents have rather similar shapes in the energy inter- 
val 0.1 GeV < e-j, < 10 TeV due to electron synchrotron 
losses. As a consequence the best observational possibil- 
ity to discriminate between leptonic and hadronic con- 
tributions at high energies 3> 1 GeV is to measure the 
7-ray spectrum at energies higher than 10 TeV, where 
we expect these two spectra to be essentially different. 
Therefore the clear detection of a substantial flux at en- 
ergies e 7 <; 10 TeV would provide direct evidence for its 
hadronic origin. 



4. Inefficient proton acceleration model 

The set of parameters which was discussed until now is 
limited by the existing measurements. There is almost no 
freedom for a substantial change of any parameter value 
without a loss of agreement between theory and experi- 
ment. For example, the proton injection rate which deter- 
mines the number of accelerated protons and the 7r°-decay 
7-ray emission, was taken to produce the shock modifica- 
tion which in turn is required to produce a steep electron 



spectrum N c cx p~ 21 needed to fit the observed radio 
emission of SN 1006. For this injection rate 77 = 2 x 10~ 4 a 
relatively high upstream magnetic field is unavoidably re- 
quired in order to produce the observed synchrotron emis- 
sion. 

On the other hand one may dismiss the relevance of a 
deviation of the observed radio spectrum from the form 
S v cx i/ -0 - 5 , and of a likely thermal contribution to the soft 
X-ray flux around 10 17 Hz, and try to reproduce all the 
observed emissions by effects of electrons alone (e.g. Pohl 
1996; Mastichiadis & de J ager 1996 Yoshida & Yanagita 



1997 ; Aharonian & Atoyan 1999|) . The number of injected 



and accelerated protons and all the effects which they can 
produce in SN 1006 is suggested to be negligibly small in 
this extreme case. For example, the proton contribution 
to the TeV 7-ray emission should be at least an order of 
magnitude lower compared with the previous case. At such 
low proton numbers the nonlinear shock modification due 
to their backreaction is negligible, and shock acceleration 
takes place in the test particle limit. 

In order to make clear how different the required set 
of relevant physical parameters compared to the case of 
efficient CR nucleon acceleration is, we performed a cal- 
culation which corresponds to inefficient CR acceleration. 
We use a proton injection rate r\ = 10~ 5 which yields an 
upper limit for the proton acceleration efficiency that is 
consistent with their low contribution to the 7-ray pro- 
duction. This low CR nucleon production efficiency also 
implies that the effect of magnetic field amplification due 
to CR streaming is small. Therefore we have used a typi- 
cal upstream ISM magnetic field strength Bq = 4 /j,G. We 
note that such a low value was even assumed to be required 
in the compressed downstream region - B± = 4 fiG - to 



fit the experimental data (Tanimori et al. |2001| ). All other 
ISM and SN parameters are the same as in the previous 
case of efficient CR acceleration. Since during the early 
evolutionary phase the global SNR dynamics is rather in- 
sensitive to the number of accelerated CRs, the observed 
SNR size and its expansion rate are equally well repro- 
duced in this case. 

In this test particle regime the electron spectrum, pro- 
duced by a strong unmodifi ed sh ock, can be analytically 
approximated by (Berezhko 1996| ) 



N c ctp 7 exp[-K(p)/«;(p^ ax )], 



(19) 



with 7 ss 2, where p^nax i s the maximum CR momentum, 
which is determined by geometrical factors. Any essential 
deviation from 7 rj 2 would be inconsistent with diffu- 
sive shock acceleration which is believed to be the main 
process of energetic particle production in SNRs. For a 
Bohm-type diffusion coefficient, sap, the cutoff region 
of this spectrum has the same form N c cx exp(— p/p^ 



as it was suggested by Mast ichiadis & de Jager (1996) 
and by Tanimori et al. ( 2001 ). At the same time, cf. Eq. 
( |l9|) , the spectrum of the electrons accelerated at a strong 
unmodified shock is essentially harder compared to what 
was used (7 = 2. 2) to fit the radio, X-ray and 7-ray data 
(Tanimori et al. 2001 ). It is clear that the spectrum cf. 
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Fig. 6. Observed total radio flux of SN 1006 (Reynolds 



1996) as a function of frequency with model spectra su- 



perimposed. Solid and dashed curves correspond in this 
figure to the high and low proton injection/acceleration 
efficiency, respectively. 




Fig. 7. Radial dependence of the 7-ray brightness for the 
7-ray energy e 7 = 3 TeV. Thick and thin curves corre- 
spond to the high and low proton injection/acceleration 
efficiency, respectively. 



Eq. (|l^) , with 7 = 2 which has the same number of GeV 
electrons to produce the same radio flux, has an order of 
magnitude more electrons with e c = 10 TeV compared 
with the case 7 = 2.2. Therefore, to fit the radio, X-ray, 
and 7-ray data with this harder electron spectrum, we 
need a large downstream magnetic field B± ps 13 [iG. 
Since even for 77 = 10~ J nucleon acceleration can pro- 
duce a strong Alfvenic wave field, the upper limit for the 
downstream field strength equals 16 /iG. To this ex- 
tent B± w 13 fj,G can still be consistent wit h a 4 /iG 
ISM magnetic field as used by Tanimori et al.( |200l| ), al- 
though these authors used 4 /iG also downstream. Since 
the energies e e <; 10 13 eV of electrons, whose synchrotron 
emission in the given magnetic field corresponds to X-ray 
energies e„ > 5 keV, should be in the exponential cutoff 
region, the required value of the maximum CR momen- 
tum is Pmax = 10 4 m p c. If we assume the same momentum 
dependence k oc p, this maximum CR momentum is con- 
sistent with the diffusion coefficient 



K(p) = 5.3fv B ohm(p), 



(20) 



which exceeds the Bohm limit by a factor 5.3 for an up- 
stream field of 4 ^G. 

The results for this inefficient proton injec- 
tion/acceleration case are also presented in Fig.^-^. 
One can see from Fig.|| that, due to the lower magnetic 
field value, one needs about five times more accelerated 
electrons compared with the previous case to fit the radio 
and X-ray data (see Fig.||). 

In the inefficient case the 7-ray production is dom- 
inated by the electron contribution (Fig.Q). Since syn- 
chrotron losses are not important, the predicted IC spec- 
trum is essentially harder for 7-ray energies between about 
10 9 and 10 11 eV than in the efficient proton acceleration 
model. The differential 7-ray CANGAROO spectrum ap- 
pears in better agreement with the 7r°-decay emission of 



the efficient model than with the inefficient model IC pre- 
diction (Fig.|). 

In Fig.|6| the calculated synchrotron fluxes are com- 
pared with the experimental data in the radio range. One 
can see that the steeper spectrum S v oc z/~ - 54 which cor- 
responds to the efficient CR acceleration case gives a bet- 
ter fit to the radio data than the spectrum S v oc v~ - 5 cor- 
responding to the case of inefficient proton injection, even 
though the experimental accuracy does not very clearly 
distinguish between these to variants. Note that a = 0.54 
is the average value of the power-law index within the 
frequency range shown in FigJ|. In fact, due to the con- 
cave shape of the electron spectrum, the index slightly 
decreases with increasing frequency, with a = 0.56 and 
a — 0.51 at the lowest and largest frequency, respectively. 

A similar situation is encountered in the X-ray range. 
An essentially steeper electron spectrum in the case of 
efficient CR accelerations much better fits the nonthermal 
X-ray data by lying below the total X-ray flux near the 
frequency v — 10 17 Hz which corresponds to the photon 
energy e„ = 0.4 keV (FigJ|). 

At given number of protons and magnetic field 
strength an electron to proton ratio K op = 0.04 is required 
in the inefficient case, in order to reproduce the observed 
emission. 

Nucleonic CRs absorb in this case only 



E c = 0.008£ sn = 2.4 x 10 49 erg, 



(21) 



which is more than an order of magnitude less than the 
canonically required 0.1£' sn for the Galactic CRs. The 
electron CR component contains E° w 10 48 erg. This is a 
f actor 7 less than what has been obtained by Dyer et al. 
( 2001 ) corresponding to their lower downstream magnetic 
field. 
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It is important to note that the 7-ray brightness of the 
remnant 



(22) 



it was already pointed out (Aharonian 1999; Aharonian & 
Atoyan |l999| Berezhko et al. |1999|; Berezhko et al. pOOl, 



J 7 (e 7 ,r) oc / dxq 7 (ej,r,x) 



has a very different dependence upon the distance from 
its center r in the two considered cases. The integration 
of the 7-ray emissivity g 7 (which is the number of pho- 
tons emitted per unit volume per second in the range e 7 
to e 7 + de-y) in the above expression is performed along 
the line of sight which intersects the visible remnant sur- 
face at the projected distance r from its center. In the 
case of inefficient proton injection TeV 7-rays are pro- 
duced by the electrons with the highest energies of their 
spectrum e > e^ ax . These electrons almost uniformly oc- 
cupy the entire volume between the shock surface R s and 
the contact discontinuity R p (which separates the ejecta 
and swept-up matter) with only a slight increase of their 
density toward the shock position. Depending on the tur- 
bulence level near the contact discontinuity they also par- 
tially fill the ejecta volume r < R p . Since the upstream 
distribution of these electrons is characterized by a dif- 
fusive length which is of the order of 0.1i? s , the region 
r > R s provides a noticeable contribution to the 7-ray 
brightness. Therefore the TeV 7-ray brightness J 7 (r) will 
have a peak value in the region R p < r < R s if particle 
penetration through the contact discontinuity is strongly 
suppressed, and at the center of the remnant r = in the 
opposite case. A significantly different radial profile J 7 (r) 
is expected in the case of efficient proton production. TeV 
7-rays are produced in this case by protons whose energy is 
considerably lower than the cutoff energy e max . Therefore 
they are concentrated within the thin region of thickness 
Ar w R s /(3a) just behind the shock. The swept-up gas 
has a similar distribution. Since the 7r°-decay 7-ray lumi- 
nosity is proportional to the product of the CR density 
and the gas number density, the 7-ray brightness J 7 (r) 
is expected to have a sharp peak close to the shock edge 
r ~ R s with a considerable decrease at small distances 
r <C -R s - The morphological effect of renormalization is to 
limit this emission to two polar caps where the interstellar 
magnetic field is parallel to the shock normal. 

The 7-ray brightness calculated for the two different 
cases is presented in Fig.^for the 7-ray energy e 7 = 3 TeV. 
One can see the essentially different radial dependence 
J 7 (r) for the two cases. The 7-rays of hadronic origin 
are expected to be concentrated in the postshock region, 
whereas the IC 7-ray brightness, calculated for quite a 
large level of turbulence near the contact discontinuity 
(see Berezhko & Volk 2000 for details) , is almost uniformly 
distributed across the visible disk of the remnant. As it is 
seen from Fig. 7, the IC high energy emission of the ineffi- 
cient case produces a noticeable halo around the shock - 
for SN1006 its size is about 0.1 degrees - that can in prin- 
ciple be used to discriminate between the 7-ray emissions 
of hadronic and leptonic origin. 

The reported 7-ray flux was detected from the same 
outer part of SN 1006 which shows the radio-emission. As 



9; Berezhko et al. li 
this can be considered as an observational argument fa- 
voring a strong role for the nuclear CR component. 



5. Summary and Conclusions 

The nonlinear kinetic model for CR acceleration in SNRs 
has been applied to SN 1006 in order to explain its ob- 
served properties. We have used stellar ejecta parameters 
M ej = 1.4M , k = 7, distance d = 1.8 kpc, and ISM 
number density iVn = 0.3 cm~ 3 from X-ray and optical 
imaginary of SN 1006. 

For these parameters an explosion energy E sn = 3 x 
10 erg is required to fit the observed size R s and expan- 
sion speed V s which are determined by the ratio E sn /N-n. 

The number of accelerated electrons required to fit the 
radio and X-ray emission of SN 1006, and correspondingly 
the role of accelerated protons in 7-ray production, de- 
pends essentially on the magnetic field value Bq. 

It was demonstrated that for low magnetic field B = 
4 fiG all the observed emissions can be dominated by the 
electron contribution. Protons are then assumed to be in- 
jected into the acceleration much less efficiently than elec- 
trons. For this test particle case the lowest permitted value 
of the electron to proton ratio is K cp = 0.04. It exceeds the 
canonical value 0.01 observed in situ in the neighborhood 
of the Solar System for the Galactic CRs. The maximum 
energy of accelerated CRs and their total energy content 
in this case are only e max ~ 10 13 eV and E c < 3 x 10 49 erg 
respectively. These numbers are too low for such SNRs to 
be considered as the main sources of the nucleonic Galactic 
CRs. 

If CRs in SN 1006 are produced due to the diffusive 
shock acceleration process, then even in the case of inef- 
ficient proton injection quite a large but plausible down- 
stream magnetic field B^ rs 13 /iG is required to fit the 
data. It is several times larger than assumed in a simple 



estimate by Tanimori et al. (2001), because the shock pro- 



duces in this case an electron spectrum N c oc e~ 2 which 
is significantly harder than the spectrum assumed in that 
estimate. 

The existing SNR data are better approximated if a 
significantly larger upstream magnetic field value Bq = 
20 /iG and a physically much more plausible, efficient nu- 
cleoli injection rate are assumed. Such an ion injection rate 
is estimated from injection theory, and consistent with the 
observed radio spectral index. The required magnetic field 
strength, that is significantly higher than the rms value 
5 /iG in the ISM, might be the result of non-linear ampli- 
fication near the SN shock by the CR acceleration process 
itself. 

The results for this case of efficient proton acceleration 
are then as follows: 

Wc find that after adjustment of the predictions of the 
nonlinear spherically-symmetric model by a renormaliza- 
tion of the number of accelerated nuclear CRs to take 
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account of the large area of quasiperpendicular shock re- 
gions in a SNR, good consistency with all observational 
data can be achieved, including the reported TeV 7-ray 
flux. The 7r°-decay 7-ray flux produced by the nuclear CR 
component exceeds the flux of 1C 7-rays generated by the 
electronic CR component at all energies above about 100 
MeV. The theory did not make use of any knowledge to be 
derived from 7-ray measurements. Therefore the reported 
TeV flux from SN 1006 supports the idea that the nuclear 
CR component is indeed produced in SNRs. The 7r°-decay 
7-ray flux comes from two polar caps of the remnant. 

The maximum energy of accelerated protons e max = 
3 x 10 14 eV and their total energy content E c w 3 x 
10 50 erg, reproduced in this case, are consistent with the 
requirements for the Galactic CR sources. The electron to 
proton ratio of K cp = 1.5 x 10~ 3 , on the other hand, is 
lower than the canonical value 0.01. 

Comparing the case of efficient proton acceleration 
with the inefficient proton acceleration case, we see that 
the expected 7r°-decay 7-ray flux oc e" 1 extends up 
to almost 100 TeV, whereas the IC 7-ray flux reaches less 
than about 10 TeV. Therefore the detection of 7-ray emis- 
sion above 10 TeV would imply evidence for a hadronic 
origin. 

We therefore conclude that the analysis of SN 1006 on 
the basis of overall SNR dynamics and nonlinear diffusive 
shock acceleration theory results in a picture where the 
nuclear component is strongly accelerated, consistent with 
all data for this SNR. 

The ratio of accelerated electrons to protons comes out 
lower by a factor of about 6 than the canonical ratio of 
0.01, generally deduced from observations of the Galactic 
CRs. If SN 1006 is typical for the nuclear Galactic CR 
source population, then other sources like young Pulsars 
or Pulsar Nebulae must also significantly contribute to 
the Galactic CR electron population. The Crab Nebula is 
a point in case. Since the majority of Galactic Supernovae 
is not of the type la to which SN 1006 belongs, but rather 
occurs as a consequence of core collapse events, there is 
the alternative possibility that core collapse Supernovae 
produce a higher electron to proton ratio. 

Phenomenological studies of the 7-ray emission from 
SN 1006 on the basis of the observed synchrotron emission 
have preferred a dominance of IC emission from electrons 
in the 7-ray part of the emission spectrum. This scenario 
can not explain the apparent bipolar morphology inferred 
from the existing 7-ray measurements and yields a much 
less convincing approximation to the radio and X-ray syn- 
chrotron spectrum. In addition, as we have shown, the IC 
emission should reach at best about 10 TeV. Apart from 
future detailed determinations of the 7-ray morphology 
one therefore needs to precisely measure the 7-ray flux at 
energies e 7 between 100 GeV and 100 TeV, and preferably 
even from 10 MeV upwards. The detailed spectra and in 
particular the existence of 7-rays with very high energies 
above 10 TeV should allow a confirmation, or a rejection, 
of CR nucleon production in SN 1006 with an acceleration 



efficiency that is consistent with the requirements on the 
Galactic CR energy budget. 
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